Identification of genetic variation influencing apolipoprotein E levels: Follow-up of genome-wide linkage scans in the GENOA study

Apolipoprotein E (ApoE) plays a major role in the metabolism of high density and low density lipoproteins (HDL and LDL). Its common protein isoforms (E2, E3, E4) are risk factors for coronary artery disease (CAD) and explain between 16 to 23% of the inter-individual variation in plasma apoE levels. Linkage analysis has been completed for plasma apoE levels in the GENOA study (Genetic Epidemiology Network of Atherosclerosis). After stratification of the population by lipoprotein levels and body mass index (BMI) to create more homogeneity with regard to biological context for apoE levels, Hispanic families showed significant linkage on chromosome 17q for two strata (LOD=2.93 at 104 cM for a low cholesterol group, LOD=3.04 at 111 cM for a low cholesterol, high HDLC group). Replication of 17q linkage was observed for apoB and apoE levels in the unstratified Hispanic and African-American populations, and for apoE levels in African-American families. Replication of this 17q linkage in different populations and strata provides strong support for the presence of gene(s) in this region with significant roles in the determination of inter-individual variation in plasma apoE levels. Through a positional and functional candidate gene approach, ten genes were identified in the 17q linked region, and 62 polymorphisms in these genes were genotyped in the GENOA families. Association analysis was performed with FBAT, GEE, and variance-component based tests followed by conditional linkage analysis. Association studies with partial coverage of TagSNPs in the gene coding for apolipoprotein H (APOH) were performed, and significant results were found for 2 SNPs (APOH_20951 and APOH_05407) in the Hispanic low cholesterol strata accounting for 3.49% of the inter-individual variation in plasma apoE levels. Among the other candidate genes, we identified a haplotype block in the ACE1 gene that contains two major haplotypes associated with apoE levels as well as total cholesterol, apoB and LDLC levels in the unstratified Hispanic population. Identifying genes responsible for the remaining 60% of inter-individual variation in plasma apoE level, will yield new insights into the understanding of genetic interactions involved in the lipid metabolism, and a more precise understanding of the risk factors leading to CAD. ^

The recombinant adeno-associated virus vector (rAAV2)-mediated apolipoprotein B mRNA-specific hammerhead ribozyme: a self-complementary AAV2 vector improves the gene expression.

BACKGROUND: In humans, overproduction of apolipoprotein B (apoB) is positively associated with premature coronary artery diseases. To reduce the levels of apoB mRNA, we have designed an apoB mRNA-specific hammerhead ribozyme targeted at nucleotide sequences GUA6679 (RB15) mediated by adenovirus, which efficiently cleaves and decreases apoB mRNA by 80% in mouse liver and attenuates the hyperlipidemic condition. In the current study, we used an adeno-associated virus vector, serotype 2 (AAV2) and a self-complementary AAV2 vector (scAAV2) to demonstrate the effect of long-term tissue-specific gene expression of RB15 on the regulation apoB mRNA in vivo. METHODS: We constructed a hammerhead ribozyme RB15 driven by a liver-specific transthyretin (TTR) promoter using an AAV2 vector (rAAV2-TTR-RB15). HepG2 cells and hyperlipidemic mice deficient in both the low density lipoprotein receptor and the apoB mRNA editing enzyme genes (LDLR-/-Apobec1-/-; LDb) were transduced with rAAV2-TTR-RB15 and a control vector rAAV-TTR-RB15-mutant (inactive ribozyme). The effects of ribozyme RB15 on apoB metabolism and atherosclerosis development were determined in LDb mice at 5-month after transduction. A self-complementary AAV2 vector expressing ribozyme RB15 (scAAV2-TTR-RB15) was also engineered and used to transduce HepG2 cells. Studies were designed to compare the gene expression efficiency between rAAV2-TTR-RB15 and scAAV2-TTR-RB15. RESULTS: The effect of ribozyme RB15 RNA on reducing apoB mRNA levels in HepG2 cells was observed only on day-7 after rAAV2-TTR-RB15 transduction. And, at 5-month after rAAV2-TTR-RB15 treatment, the apoB mRNA levels in LDb mice were significantly decreased by 43%, compared to LDb mice treated with control vector rAAV2-TTR-RB15-mutant. Moreover, both the rAAV2-TTR-RB15 viral DNA and ribozyme RB15 RNA were still detectable in mice livers at 5-month after treatment. However, this rAAV2-TTR-RB15 vector mediated a prolonged but low level of ribozyme RB15 gene expression in the mice livers, which did not produce the therapeutic effects on alteration the lipid levels or the inhibition of atherosclerosis development. In contrast, the ribozyme RB15 RNA mediated by scAAV2-TTR-RB15 vector was expressed immediately at day-1 after transduction in HepG2 cells. The apoB mRNA levels were decreased 47% (p = 0.001), compared to the control vector scAAV2-TTR-RB15-mutant. CONCLUSION: This study provided evidence that the rAAV2 single-strand vector mediated a prolonged but not efficient transduction in mouse liver. However, the scAAV2 double-strand vector mediated a rapid and efficient gene expression in liver cells. This strategy using scAAV2 vectors represents a better approach to express small molecules such as ribozyme.

Origins and evolution of VNTR loci: The apolipoprotein B 3$\sp\prime$ VNTR

I studied the apolipoprotein (apo) B 3$\sp\prime$ variable number tandem repeat (VNTR) and did computer simulations of the stepwise mutation model to address four questions: (1) How did the apo B VNTR originate? (2) What is the mutational mechanism of repeat number change at the apo B VNTR? (3) To what extent are population and molecular level events responsible for the determination of the contemporary apo B allele frequency distribution? (4) Can VNTR allele frequency distributions be explained by a simple and conservative mutation-drift model? I used three general approaches to address these questions: (1) I characterized the apo B VNTR region in non-human primate species; (2) I constructed haplotypes of polymorphic markers flanking the apo B VNTR in a sample of individuals from Lorrain, France and studied the associations between the flanking-marker haplotypes and apo B VNTR size; (3) I did computer simulations of the one-step stepwise mutation model and compared the results to real data in terms of four allele frequency distribution characteristics.^ The results of this work have allowed me to conclude that the apo B VNTR originated after an initial duplication of a sequence which is still present as a single copy sequence in New World monkey species. I conclude that this locus did not originate by the transposition of an array of repeats from somewhere else in the genome. It is unlikely that recombination is the primary mutational mechanism. Furthermore, the clustered nature of these associations implicates a stepwise mutational mechanism. From the high frequencies of certain haplotype-allele size combinations, it is evident that population level events have also been important in the determination of the apo B VNTR allele frequency distribution. Results from computer simulations of the one-step stepwise mutation model have allowed me to conclude that bimodal and multimodal allele frequency distributions are not unexpected at loci evolving via stepwise mutation mechanisms. Short tandem repeat loci fit the stepwise mutation model best, followed by microsatellite loci. I therefore conclude that there are differences in the mutational mechanisms of VNTR loci as classed by repeat unit size. (Abstract shortened by UMI.) ^

A cladistic analysis of apolipoprotein B gene variation and plasma lipid and apolipoprotein B levels, using familial data

Any functionally important mutation is embedded in an evolutionary matrix of other mutations. Cladistic analysis, based on this, is a method of investigating gene effects using a haplotype phylogeny to define a set of tests which localize causal mutations to branches of the phylogeny. Previous implementations of cladistic analysis have not addressed the issue of analyzing data from related individuals, though in human studies, family data are usually needed to obtain unambiguous haplotypes. In this study, a method of cladistic analysis is described in which haplotype effects are parameterized in a linear model which accounts for familial correlations. The method was used to study the effect of apolipoprotein (Apo) B gene variation on total-, LDL-, and HDL-cholesterol, triglyceride, and Apo B levels in 121 French families. Five polymorphisms defined Apo B haplotypes: the signal peptide Insertion/deletion, Bsp 1286I, XbaI, MspI, and EcoRI. Eleven haplotypes were found, and a haplotype phylogeny was constructed and used to define a set of tests of haplotype effects on lipid and apo B levels.^ This new method of cladistic analysis, the parametric method, found significant effects for single haplotypes for all variables. For HDL-cholesterol, 3 clusters of evolutionarily-related haplotypes affecting levels were found. Haplotype effects accounted for about 10% of the genetic variance of triglyceride and HDL-cholesterol levels. The results of the parametric method were compared to those of a method of cladistic analysis based on permutational testing. The permutational method detected fewer haplotype effects, even when modified to account for correlations within families. Simulation studies exploring these differences found evidence of systematic errors in the permutational method due to the process by which haplotype groups were selected for testing.^ The applicability of cladistic analysis to human data was shown. The parametric method is suggested as an improvement over the permutational method. This study has identified candidate haplotypes for sequence comparisons in order to locate the functional mutations in the Apo B gene which may influence plasma lipid levels. ^

THE NEW ROLE OF PROPROTEIN CONVERTASE SUBTILISIN/KEXIN TYPE 9: A CONNECTION OF PROPROTEIN CONVERTASE SUBTILISIN/KEXIN TYPE 9, APOLIPOPROTEIN B, and AUTOPHAGY

Plasma low-density lipoprotein (LDL) levels are positively correlated with the incidence of coronary artery disease. In the circulation, the plasma LDL clearance is mainly achieved by the uptake via LDL receptor (LDLR). Proprotein convertase subtilisin/kexin type 9 (PCSK9) is a newly discovered gene, playing an important role in LDL metabolism. Gain-of-function mutations of PCSK9 lead to hypercholesterolemia and loss-of-function mutations of PCSK9 are associated with decrease of LDL cholesterol. The effects of PCSK9 on cholesterol levels are the consequence of a strong interaction between the catalytic domain of PCSK9 and epidermal growth factor-like repeat A (EGF-A) domain of LDLR on the cell surface of hepatocytes. This PCSK9/LDLR complex enters the cell via endocytosis, where both PCSK9 and LDLR are removed via the lysosome pathway, resulting in decreased levels of LDLR and accumulation of LDL in the plasma. However, whether this is the exclusive function of PCSK9 on LDL metabolism was challenged by us; we observed PCSK9 interacted with apolipoprotein B (apoB) and increased apoB production, irrespective of the LDLR. ApoB is the primary structure protein of LDL particle and it also serves as the ligand for the LDL receptor. There is ample evidence showing that the levels of apoB are a better indicator for heart disease than either total cholesterol or LDL cholesterol levels. We used a second-generation adenoviral vector to overexpress PCSK9 (Ad-PCSK9) in wild-type C57BL/6 and LDLR deficient mice (Ldlr-/- and Ldlr-/-Apobec1-/-). Our study revealed that overexpression of PCSK9 promoted the production and secretion of apoB in the form of very-low density lipoprotein (VLDL), which is the precursor of LDL, in the 3 mouse models studied (C57BL/6J, Ldlr-/-, and Ldlr-/-Apobec1-/-). The increased apoB production in mice was regulated at post-transcriptional levels, since there was no difference in apoB mRNA levels between mice treated with Ad-PCSK9 and control vector Ad-Null. By using pulse-chase experiment on primary hepatocytes, we showed that overexpression of PCSK9 increased the secretion of apoB, independent of LDLR. In the circulation, we showed that PCSK9 was associated with LDL particles. By using 3 different protein–protein interaction assays of co-immunoprecipitation, mammalian two-hybrid system, and in situ proximity ligation assay, we demonstrated a direct protein–protein interaction between PCSK9 and apoB. The impact of this interaction inhibited the physiological removal process of apoB via autophagosome/lysosome pathway in an LDLR-independent fashion, resulting in increased production and secretion of apoB-containing lipoproteins. The significance of this process was shown in the Pcsk9 knockout mice in the background of Ldlr-/-Apobec1-/- mice (triple knockout mice); in the absence of Pcsk9 (triple knockout mice) the levels of cholesterol, triacylglycerol, and apoB decreased significantly in comparison to that of Ldlr-/-Apobec1-/- mice. Taken together, our study demonstrated a direct intracellular interaction of PCSK9 with apoB, resulting in the inhibition of apoB degradation via the autophagosome/lysosome pathway independent of LDLR. This discovery provides a new concept of the importance of PCSK9 and suggests new approaches for the therapeutic intervention of hyperlipidemia.

NOVEL THERAPEUTIC TARGETS IDENTIFIED BY HIGH-THROUGHPUT TECHNOLOGIES FOR TRIPLE-NEGATIVE BREAST CANCER

Triple-negative breast cancers (TNBC) are characterized by the lack of or reduced expression of the estrogen and progesterone receptors, and normal expression of the human epidermal growth factor receptor 2. The lack of a well-characterized target for treatment leaves only systemic chemotherapy as the mainstay of treatment. Approximately 60-70% of patients are chemosensitive, while the remaining majority does not respond. Targeted therapies that take advantage of the unique molecular perturbations found in triple-negative breast cancer are needed. The genes that are frequently amplified or overexpressed represent potential therapeutic targets for triple-negative breast cancer. The purpose of this study was to identify and validate novel therapeutic targets for triple-negative breast cancers. 681 genes showed consistent and highly significant overexpression in TNBC compared to receptor-positive cancers in 2 data sets. For two genes, 3 of the 4 siRNAs showed preferential growth inhibition in TNBC cells. These two genes were the low density lipoprotein receptor-related protein 8 (LRP8) and very low-density lipoprotein receptor (VLDLR). Exposure to their cognate ligands, reelin and apolipoprotein E isoform 4 (ApoE4), stimulated the growth of TNBC cells in vitro. Suppression of the expression of either LRP8 or VLDLR or exposure to RAP (an inhibitor of ligand binding to LRP8 and VLDLR) abolished this ligand-induced proliferation. High-throughput protein and metabolic arrays revealed that ApoE4 stimulation rescued TNBC cells from serum-starvation induced up-regulation of genes involved in lipid biosynthesis, increased protein expression of oncogenes involved in the MAPK/ERK and DNA repair pathways, and reduced the serum-starvation induction of biochemicals involved in oxidative stress response and glycolytic metabolism. shLRP8 MDA-MB-231 xenografts had reduced tumor volume, in comparison to parental and shCON xenografts. These results indicate that LRP8-APOE signaling confers survival advantages to TNBC tumors under reduced nutrient conditions and during cellular environmental stress. We revealed that the LRP8-APOE receptor-ligand system is overexpressed in human TNBC. We also demonstrated that this receptor system mediates a strong growth promoting and survival function in TNBC cells in vitro and helps to sustain the growth of MDA-MD-231 xenografts. We propose that inhibitors of LRP8-APOE signaling may be clinically useful therapeutic agents for triple-negative breast cancer.

Gene by BMI Interactions Influencing C-Reactive Protein Levels in European-Americans

C-Reactive Protein (CRP) is a biomarker indicating tissue damage, inflammation, and infection. High-sensitivity CRP (hsCRP) is an emerging biomarker often used to estimate an individual’s risk for future coronary heart disease (CHD). hsCRP levels falling below 1.00 mg/l indicate a low risk for developing CHD, levels ranging between 1.00 mg/l and 3.00 mg/l indicate an elevated risk, and levels exceeding 3.00 mg/l indicate high risk. Multiple Genome-Wide Association Studies (GWAS) have identified a number of genetic polymorphisms which influence CRP levels. SNPs implicated in such studies have been found in or near genes of interest including: CRP, APOE, APOC, IL-6, HNF1A, LEPR, and GCKR. A strong positive correlation has also been found to exist between CRP levels and BMI, a known risk factor for CHD and a state of chronic inflammation. We conducted a series of analyses designed to identify loci which interact with BMI to influence CRP levels in a subsample of European-Americans in the ARIC cohort. In a stratified GWA analysis, 15 genetic regions were identified as having significantly (p-value < 2.00*10-3) distinct effects on hsCRP levels between the two obesity strata: lean (18.50 kg/m2 < BMI < 24.99 kg/m2) and obese (BMI ≥ 30.00 kg/m2). A GWA analysis performed on all individuals combined (i.e. not a priori stratified for obesity status) with the inclusion of an additional parameter for BMI by gene interaction, identified 11 regions which interact with BMI to influence hsCRP levels. Two regions containing the genes GJA5 and GJA8 (on chromosome 1) and FBXO11 (on chromosome 2) were identified in both methods of analysis suggesting that these genes possibly interact with BMI to influence hsCRP levels. We speculate that atrial fibrillation (AF), age-related cataracts and the TGF-β pathway may be the biological processes influenced by the interaction of GJA5, GJA8 and FBXO11, respectively, with BMI to cause changes in hsCRP levels. Future studies should focus on the influence of gene x bmi interaction on AF, age-related cataracts and TGF-β.

GENETIC EVENTS PREDISPOSING TO WILMS' TUMOR

Wilms' tumor (WT) is a childhood embryonic tumor of the kidney. In some cases, WT has been associated with a chromosome deletion in the region 11p13. The majority of WT cases, however, have normal karyotypes with no discernable deletions or rearrangements of chromosome 11.^ To study the genetic events predisposing to the development of WT, I have used a number of gene markers specific for chromosome 11. Gene probes for human catalase and apolipoprotein A1 were localized to chromosome 11 by in situ hybridization. A number of other probes previously mapped to chromosome 11 were also used. Nine WT patients who were heterozygous for at least one 11p marker were shown to lose heterozygosity in their tumor DNA. Gene dosage experiments demonstrated that two chromosomes 11 were present although loss of heterozygosity had occurred in all but two cases. By using gene probes from the short and long arms of chromosome 11, I discerned that loss of heterozygosity was due to somatic recombination in four cases, chromosome deletion in two cases, and chromosome loss and reduplication or somatic recombination in these cases. Examination of DNAs from the parents of six of these patients indicated that the alleles that were lost in tumor tissues were alleles inherited from the mother. In sporadic WT cases one would expect the loss of alleles to be random. These data suggest that the loss of alleles resulting in the development of WT is not a random event, however, the significance of this is not known. ^

An ordinal logistic regression model with misclassification of the outcome variable and categorical covariate

Ordinal outcomes are frequently employed in diagnosis and clinical trials. Clinical trials of Alzheimer's disease (AD) treatments are a case in point using the status of mild, moderate or severe disease as outcome measures. As in many other outcome oriented studies, the disease status may be misclassified. This study estimates the extent of misclassification in an ordinal outcome such as disease status. Also, this study estimates the extent of misclassification of a predictor variable such as genotype status. An ordinal logistic regression model is commonly used to model the relationship between disease status, the effect of treatment, and other predictive factors. A simulation study was done. First, data based on a set of hypothetical parameters and hypothetical rates of misclassification was created. Next, the maximum likelihood method was employed to generate likelihood equations accounting for misclassification. The Nelder-Mead Simplex method was used to solve for the misclassification and model parameters. Finally, this method was applied to an AD dataset to detect the amount of misclassification present. The estimates of the ordinal regression model parameters were close to the hypothetical parameters. β1 was hypothesized at 0.50 and the mean estimate was 0.488, β2 was hypothesized at 0.04 and the mean of the estimates was 0.04. Although the estimates for the rates of misclassification of X1 were not as close as β1 and β2, they validate this method. X 1 0-1 misclassification was hypothesized as 2.98% and the mean of the simulated estimates was 1.54% and, in the best case, the misclassification of k from high to medium was hypothesized at 4.87% and had a sample mean of 3.62%. In the AD dataset, the estimate for the odds ratio of X 1 of having both copies of the APOE 4 allele changed from an estimate of 1.377 to an estimate 1.418, demonstrating that the estimates of the odds ratio changed when the analysis includes adjustment for misclassification. ^